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ABSTRACT 

Numerical simulations of a planar turbulent 
wall jet and a planar VTOL upwash fountain have 
been performed. These are three-dimensional simu- 
lations and they resolve large scale unsteady 
motions in the flows. The wall jet simulation 
shows good agreement with experimental data and is 
presented to verify the simulation methodology. 
Simulation of the upwash fountain predicts ele- 
vated shear stress and a half -velocity width 
spreading rate of 33% which agrees well with ex- 
periment. Turbulence mechanisms which contribute 
to the enhanced spreading rate are examined. 

1 . INTRODUCTION 

The desire to have a VTOL supersonic fighter 
has increased in recent years as there are many 
advantages in having an aircraft that can be based 
Independent of conventional runways. VTOL capa- 
bility is provided by some combination of downward 
thrusting jets. In ground effect these jets pro- 
duce fluid dynamical problems that are not typi- 
cally encountered in conventional aircraft. A 
complete list of ground effects problems would be 
rather large and would include, for example, the 
Reynolds number scaling of the "suck -down" effect, 
the enhanced spreading rate of the upwash foun- 
tain, hot gas reingestion of the fountain or 
ground vortex fluid, and aircraft stability prob- 
lems due to interaction with the fountain or the 
ground vortex. 

The key to an understanding and a predictive 
capability in many of the above problems is in the 
ability to understand the turbulent mixing. The 
suck^down effect and the fountain's spreading rate 
are almost purely turbulence problems; the ground 
vortex also depends on inviscid and bouyancy phe- 
nomena. These are complex turbulent flows, invol- 
ving combinations of "turbulence modifiers" that 
are rarely encountered in other applications and, 
therefore, it will be difficult to develop ade- 
quate models for these flows. 

This paper focuses on the turbulence in the 
upwash fountain, sketched in Figure 1, which Is 
important for several practical reasons. The 
upwash is hot, it strikes the underside of the 
aircraft, it contributes to lift, and it may lead 
to exhaust gas reingestion. From a scientific 
standpoint the upwash spreading rate is an anomaly 
which begs to be explained. The combination of 
8 t reamline curvature and velocity gradient can 
have a pronounced effect on turbulent mixing 
(Refs. 1, 2, 3]. Present understanding indicates 
that the combination of curvature and velocity 
gradient, like that in a boundary layer on a con- 
vex surface, should stabilize the turbulence in 
most of the turning region of the fountain. How- 
ever, just the opposite occurs and the shear 
stresses and spreading rate are two to three times 


greater than in a conventional jet. There is a 
near -wall region which has shear of the opposite 
sign, and for which the curvature -velocity gra- 
dient combination effect should be destabi- 
lizing. This region is small in the wall jets and 
in the fountain, but in the collision zone it 
encompasses much of the reversed-f low zone under 
the fountain. However, it seems unlikely that 
this relatively small region of the flow could 
energize the turbulence in all of the fountain. 

The high stresses and spreading rate are 
underpredicted by most classes of current turbu- 
lence models: the k-e, the algebraic Reynolds 

stress, and the Reynolds stress transport models 
[Ref. 4 and Launder, private communication]. 
These models seem to be insensitive to the turbu- 
lence mechanisms in the fountain's turning re- 
gion. For obvious reasons it is desirable to have 
a model that gives improved results for this 
flow. Additionally, It is clear that these poten- 
tially significant mechanisms are not represented 
in the modeling of other flows, where their ab- 
sence, not being as critical, has escaped detec- 
tion. The essential feature of the upwash foun- 
tain Is colliding wall -bounded shear layers with 
strong curvature at the collision point. A two- 
dimensional boundary layer with a reversed flow 
region has similar features at the separation 
point. 

Finally, it is worth noting that experimental 
studies of the fountain are difficult to per- 
form. The flow is highly unsteady with frequent 
flow direction reversal in the critical region 
where the wall jets collide. Laser instrumenta- 
tion would seem to be necessary for this work, but 
has only recently been employed [Ref. 6]. Fur- 
thermore, the large pressure gradients in a region 
of unsteadiness Indicate that pressure fluctua- 

tions may be important and these cannot be mea- 
sured at present. The pressure-strain term in the 
Reynolds stress transport equations has, there- 
fore, not been directly measured although it is 
thought to be important in curved flows. 

1.1 The Approach 

In this work two flows are simulated, the 
planar wall jet, chosen because it is a well docu- 
mented flow which can be used to validate the 
numerical method, and the planar upwash foun- 
tain. There are many advantages to studying the 

planar flow rather than the fountain resulting 

from radial wall jets, which more closely resem- 
bles the VTOL ground effect flow field. These 
include the ability to use periodic boundary con- 
ditions In one direction and shorter computing 

times for good statistics. 
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The simulations are done by specifying un- 
steady inflow conditions that approximate the 
turbulent wall jets. These flows then evolve in 
the strearawise direction before they are evalu- 
ated, in the case of the wall jet simulation, or 
collide with another wall jet, in the fountain 
simulation. Instantaneous and time -averaged data 
are obtained from the simulation. This makes the 
work similar to an experimental program except 
that the type and quantity of data available can 
be greater. To date single point correlations of 
velocities, velocity gradients, pressure and den- 
sity have been obtained. These terms give us, 
from the Reynolds —averaged standpoint, the ability 
to determine most of the mechanisms critical to 
the turbulence in the fountain. This work is 
similar to the work in Large Eddy Simulation (LES) 
reviewed by Rogallo and Moin [Ref. 7], for exam- 
ple. Many differences distinguish LES from the 
present effort, which we call Very Large Eddy 
Simulation (VLES). 

2. METHODS 

2.1 Equation Set and Numerical Algorithm 


it is assumed, are more easily modeled and less 
critical to the turbulent processes than the large 
scales. We use an eddy viscosity which is propor- 
tional to the magnitude of vorticity 


CA 


( 1 ) 


The value of C tT represents the square of 
the SGS mixing length. A constant characteristic 
grid spacing A is used, since the grid spacing in 
the x- and z -direct ions is nearly constant in the 
collision zone. The coefficient C takes the value 
of C ■ 0. 14 which is close to the values used for 
LES in References 10 and 11, except near the 
wall. At the wall C approaches zero according to 
the Von Karraan formula so that the log -law profile 
is reoroduced. 


2.2 Boundary Conditions 

There are several types of boundary condi- 
tions which must be applied in these simula- 
tions. The straight-forward ones will be dis- 
cussed first, then the more difficult inflow and 
outflow conditions will be given. 


The simulations have been performed using a 
conventional finite difference method that is 
often used for steady-state Reynolds -averaged 
calculations of aerodynamics problems. The vis- 
cous conservation equations for mass, energy, and 
momenta in three directions (commonly, the Navier- 
Stokes equations) are used for the simulations. 
They are written in a conservation law form common 
for high speed aerodynamics. The solution algo- 
rithm is the 1981 implicit -explicit MacCormack 
method [Ref. 8]. The implementation of the algo- 
rithm includes modifications to improve its effi- 
ciency for this application; these are discussed 

in Childs and Nixon [Ref. 5]. The two significant 
modifications are: 1) viscous diffusion terms are 

first order in time, as opposed to the second 
order accuracy of the original algorithm (not a 
significant disadvantage since the sub -grid -scale 
turbulence model, which dominates viscous diffu- 
sion, has no formal accuracy at all, and convec- 
tive terras are still second order in space and 

time); and 2) the time step is chosen such that 

the algorithm is explicit in the two directions 
parallel to the wall. The implicit step is only 

used normal to the wall. Due to the grids which 
are used the algorithm is entirely explicit away 
from the grid clustering at the wall* This mini- 
mizes the dissipative effects of an implicit 

method, but does not overly restrict the time 
step. 

The difficulty at solid boundaries with the 
1981 MacCormack method [Ref. 9] has not been a 
problem in the present work since the time steps 
needed to resolve the turbulence are relatively 
small and give a maximum CFL number of 20, typi- 
cally. 

The simulations are performed on rather 
coarse grids because of computer limitations and 
only the largest scales of turbulence can be simu- 
lated. Therefore, the equations are Reynolds- 
averaged and a sub-grid -scale (SGS) model is used 
for the unresolvable scales of turbulence which, 


In the z-direction the flow is periodic. 
This condition is imposed by over-writing the 
boundary points with data from the first interior 
point at the opposite side of the grid. At the 
solid wall the no-slip condition is applied to 
velocities and a zero-gradient condition is ap- 
plied to the density and energy. The wall falls 
between two grid points and these conditions imply 
a zero-gradient on temperature and pressure. 

At the top of the domain there is outflow in 
the fountain and slow flow, either in or out, on 
either side of the fountain. This boundary is 
treated with a zero-gradient extrapolation condi- 
tion on all variables. There is no effort to 
obtain meaningful data at this boundary. The 
primary concern is that there be no reflection of 
pressure waves. From analysis of outflow boundary 
conditions (c.f. Ref. 12) the imposition of pres- 
sure and extrapolation of density and velocities 
is correct for steady calculations. However, 
imposing a pressure at the top would produce pres- 
sure reflections. 

The side boundaries above the wall jets oust 
provide a small amount of inflow for entrainment 
into the jet. It would be incorrect to specify 
the inflow because that could amount to specifying 
the entrainment into the jet. A zero-gradient 
condition on the velocities permits the inflow to 
adjust as necessary to satisfy the entrainment. 
The density is set to ambient. At some point the 
pressure must be tied to ambient conditions and 
there is no other suitable place to do this since 
all other boundaries will experience significant 
pressure fluctuations. A H loosely tied pressure'^ 
condition is used, which provides a small "puU" 
towards ambient pressure, but permits the pressure 
to deviate if the interior flow so requires. The 
loosely tied pressure is computed as a weighted 
average of the local interior pressure and the 
fixed ambient pressure. Typically the weighting 
is about 30% on the ambient pressure and 70% on 
the interior pressure. A similar type of pressure 
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boundary condition was used by Boris [Ref. 13] in 
the simulation of turbulent outflow. 

The inflow conditions are clearly the most 
difficult. They must approximate the mean and 
turbulent flow in a fully developed turbulent wall 
jet. There were two guiding principles used in 
determining the inflow. Firstly, it should give 
correct values for the mean flow and basic statis- 
tics, the normal and shear stresses. Secondly, 
the unsteady inflow should "look right" when com- 
pared to flow visualization pictures, A third 
condition, on the turbulent energy spectrum, may 
be examined in the future. 

The unsteady inflow profiles are constructed 
from a combination of Chebyshev modes normal to 
the wall and Fourier modes in the z-direction and 
time. This is added to a mean inflow profile 
determined from experimental data [Ref. 14, p, 
434]. The following expression gives the unsteady 
8 treamwl 8 e velocity perturbation 

u '(y,z,t ) - U max a t>B T m (y) «in( V ± c^t) 

( 2 ) 

The Chebyshev polynomials are represented by 
^ are wavelengths, the are wave 
speeds and the are weighting coefficients. 

The lateral and vertical components are specified 
as functions of the streamwise fluctuations. 

v " - u" [± 0 X + 0 2 sin<w z) sin(t)] 

P O) 

w * u [ 0 3 sin(u) p z) sin(t)] 

The many coefficients in Equations ( 2 ) and 
(3) must satisfy certain constraints. The velo- 
city perturbation goes to zero at the wall and the 
edge of the jet, and the maximum values of the 
normal and shear stresses must agree with experi- 
ment. These conditions do not uniquely determine 
the coefficients and the condition on the "appear- 
ance" of the unsteady profile was employed to 
determine ratios between some coefficients. The 
inflow density was set to ambient and the pressure 
was extrapolated from the interior. Specification 
of the above boundary conditions completes the 
description of the method. 

2.3 Grids, Statistics and Computation Times 

Cartesian grids with stretching in the x- and 
y -directions are used. The grid is clustered at 
the wall and in the center of the domain (for 
upwash simulations) where the fountain is expected 
to form. Simulations have been run on two grids 
with different resolutions for both the wall jet 
and the fountain. The results showed grid depen- 
dence in magnitude but not in character; in par- 
ticular, the spreading rate of the fountain was 
very similar in both cases. The finer grid 

results are shown. For the fountain the grid has 
40 by 32 by 25 points in the x, y, and z direc- 
tions; for the wall jet the grid has 32 by 32 by 
32 points. 


Analysis of the turbulence is done by accumu- 
lating single point statistics of velocities, 
velocity gradients, pressure, density and 
energy. These are then processed to provide the 
turbulent correlations. The flow is two- 
dimensional in the mean and statistical averaging 
is done over time and the z-direction, which is 
the mean-invariant direction. Turbulence evalua- 
tion is based on the assumption of constant den- 
sity, which is sufficiently accurate for the pre- 
sent purposes. 

The integration times are, so far, marginally 
adequate for some correlations, which are still 
changing slowly with time. Second order correla- 
tions, which give the Reynolds stresses and pres- 
sure-strain terms, appear to be stable to within 
5% for the wall jet simulation but variations of 
20% may be seen in the fountain results. Higher 
order statistics are not sufficiently converged 
for many purposes. Variations of 50% may be seen 
before they stabilize. The reason is that higher 
order statistics are formed from small differences 
between lower order statistics. Small variations 
in the lower order statistics can produce large 
changes in the higher order correlation. The 
higher order statistics are given to show orders 
of magnitude and trends, which are firmly estab- 
lished. 

These calculations have been run on a CRAY X- 
MP. The upwash fountain was run for 90 minutes of 
CPU time and statistics were taken over the last 
45 minutes. This provided 110 units of physical 
time nondimensionalized on the acoustic speed and 
initial wall jet half-velocity thickness, or ap- 
proximately 1.7 "flow-through" time periods (time 
for the maximum-velocity fluid to transit the 
domain) • 

3. RESULTS 

The results will be presented in two 
stages. The first step is to validate the simula- 
tion, which is done for the wall Jet. The accu- 
racy of the wall j’fct simulation will indicate what 
accuracy can be expected for the fountain simula- 
tion, The second stage is to examine the results 
for the fountain and investigate the turbulence 
mechanisms. 

3.1 The Wall Jet 

The simulation is of a two-dimensional wall 
jet on a plane surface with no streamwise pressure 
gradient. The calculation was run with a mean 
maximum Inflow Mach number of 0.65, which is slow 
enough to minimize compressibility effects yet 
fast enough to maintain good computational effi- 
ciency. The wall and the "freest ream" were given 
a velocity of 10% of the inflow maximum. This 
gives clearly defined inflow and outflow boun- 
daries. Since the freestream and wall have the 
same velocity the flow has the conditions required 
for self-preservation. The added velocity is 
removed from all data analysis and the results. 
The Reynolds number based on maximum velocity, U , 
and half -velocity thickness, 20000. Tne 

data to be used for comparison were compiled by 
Launder and Rodi for the 1980-81 Stanf ord-AFOSR 
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Conference [Ref* 14] from several experiments 
referenced therein. 

For the comparison we use time-averaged 
results from the unsteady simulation* Figure 2 
gives the velocity vectors at every fourth stream- 
wise grid line. The dimensions in all two- 

dimensional figures are normalized on the half- 
velocity thickness at the inlet. The inlet velo- 
city profile is specified to match experimental 
data. Between x-1 and x*9 the profile departs 

from experiment and shows a velocity profile which 
is too steep in the outer third of the jet. By 
x*10 the velocity profile is once more in agree- 
ment with experiment, as shown in Figure 3. 
There are small errors near the wall, but this is 
where the grid resolution in the x- and z- 
directions is the poorest relative to the need, 
and most of the shear stress is carried by the SGS 
model. The growth rate of the half -velocity 
thickness (at x*10) is 0.067, which is 9% below 
the experimental value of 0.073. Although the 
mean velocity profile has stabilized at this 

point, many of the turbulent statistics have 
not. The flow is not yet self -preserving. 

In Figure 4 the resolvable components of the 
normal stresses are given. Clearly, the SGS con- 
tribution to the stresses should be positive 

(which eddy viscosity models don't guarantee) and, 
thus, the total stresses would be greater than the 
values shown. The resolvable stresses are all 
within 20 Z of the experiment and predict the 
trends that are important, including the slight 
rise in the u' and w ' near the wall. The reduc- 
tion in u' and w' very near the wall is physically 
correct behavior, although the region where this 
occurs is too thick in the simulation. This is 

likely caused by inadequate numerical resolution, 
but may also be due to the SGS model. 

The resolvable and resolvable— plus -SGS shear 
stresses are shown in Figure 5. Consider first 
only the resolvable stress. The stress in the 
outer 75% of the wall jet is well predicted, al- 
though the location and magnitude of stress at its 
maximum are slightly high. In the near-wall re- 
gion the stress is badly underpredicted. Consid- 
ering the poor resolution here this problem was 
not unexpected. When the SGS stress is included, 
the near-wall results improve and give nearly the 
correct skin friction. However, the stress in the 
outer layer is now overpredicted by about 60%. 

Figure 6 gives the tri ple ve locity correla- 
tions of u'u'u', u'v'v' and v'v'v'. These terms 
are responsible for the turbulent transport of the 
Reynolds stresses. They are important to under- 
standing the turbulence, to multi -equation models, 
and to validating the present simulations. The 
results show that, in general, the trends and 
orders of magnitude are well predicted, with t he 
exception that the near -wall region of v'v'v' is 
not. The overprediction of u'u'u' in the outer 
layer is consistent with the overpr ediction there 
0 f u V. The prediction of u'v'v' (shear stress 
transport) is the best of the three and is only in 
significant error in the near-wall region. It 
should be noted that measurement errors of triple 
products are of the order of 15% to 30%, under the 
best conditions. 


In summary for the wall jet the following are 
the major points. The resolvable normal stresses 
are predicted to within 20% of experiment, with 
the streamwise fluctuations being high and the 
vertical fluctuations being low. The shear stress 
in the outer layer is overpredicted by a signifi- 
cant amount. We do not believe this error will 

degrade the fountain simulation, for two rea- 
sons. In the fountain simulation the grid resolu- 
tion of the wall jets is less and the turbulence 
levels are lower than in the pure wall jet simula- 
tion. Elevated turbulence in the fountain is, 
more than ever, due to mechanisms in the collision 
zone. Also, the turbulence levels In the fountain 
are many times greater than In the wall jet , so 
the turbulence production In the collision zone 
over-powers the wall jet turbulence. The third 
order correlations are adequate away from the 
wall, but are underpredicted in the near-wall 
region, which is due primarily to inadequate grid 
resolution in the x- and z-directions . In gen- 
eral, the turbulence in the wall jet is suffi- 
ciently close to experiment to permit useful study 
of the upwash turbulence. 

3.2 The Upwash Fountain 

The objective of this study is to examine the 
turbulence in the upwash fountain. The simulation 
was run by colliding two turbulent wall jets and 
permitting the upwash to develop naturally. The 
wall jets are run with the same inlet conditions 
as the pure wall jet. 

For a point of comparison, the results from a 
steady state Reynolds-averaged calculation using 
a k-e model are given in Figure 7. Although not 
shown, the turbulence levels and spreading rate 
are typical for the values seen in a "normal" free 
jet, and significantly underpredict the values 
seen in upwash fountains. 

Figure 8 gives the mean velocity vectors and 
Mach number contours from the unsteady simulation, 
similar to Figure 7 from the steady calculation. 
Rapid spreading of the fountain is clearly seen. 
Initially a velocity deficit at the center of the 
fountain is seen but It quickly disappears above 
the collision zone. Streamlines, shown in Figure 
8-c, reveal a two stage collision process. First 
the wall jets separate and flow over relatively 
large reversed flow regions. Then they collide 
above the separation zone and are redirected up 
ward. The half velocity width of the fountain is 
plotted in Figure 9; it grows at a rate of 33% 
just above the collision zone. This compares well 
with experimental values in the range of 24% to 
32% [Refs. 15, 16]. The data of Gilbert [Ref. 15] 
are included and show good agreement midway up the 
fountain, but poorer agreement at the bottom and 
top. It is likely that measurement errors (hot 

wire anemometry) are significant at the bottom and 
that simulation errors are significant at the 
top. The fountain must eventually relax to a 

conventional plane jet, although this may be oc- 
curring too rapidly in the simulation. Turbulence 
near the top boundary Is damped by the dissipative 
outflow boundary conditions and coarser grid spac 
ing, which decreases the spreading rate. Large 
scale motions are also constrained by the periodi- 
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city condition because the width of the fountain 
near the top boundary is greater than the z- 
dimension of the physical domain (see Fig. 16). 

Figure 10 gives the mean velocity profiles in 
the fountain at several heights. The fountain is 
not perfectly symmetric although it should be. 
The factors which contribute to this are non- 
symmetric truncation errors in the MacCormack 
algorithm and insufficient integration time. It 
should be noted that most experimental results 
show some asymmetry and, thus, the fountain may 
easily be disturbed from symmetric. 

The streamwise normal stress is given in 
Figure 11. Data from physical experiments of 
Gilbert [Ref. 15], Kind [Ref. 16] and Saripalli 
[Ref. 6] are included. The "normal" plane free 
jet has been studied by Bradbury [Ref. 17] and his 
results are quite close to those of Gilbert. The 
simulated results are below those from Kind but 
generally above those of Gilbert. Also, a sharp 
dip in the stress is seen at the centerline in the 
present results and the results of Kind, but not 
in those of Gilbert or Saripalli. 

Contour plots of some turbulence quantities 
are given in what follows. In these, the normal- 
izing velocity scale is /2Ap, where Ap is the 
maximum pressure rise in the collision zone. This 
is approximately the maximum mean velocity of the 
colliding wall jets. The normalizing length scale 
i® half of the width of the high pressure re- 
gion. The time scale derived from the length and 
velocity scales is approximately the minimum time 
for fluid to transit through the collision zone. 
Figure 12 gives contour p lots of the normal 
stresses. The maximum for u'u" is 0.28 (note 
these are lateral fluctuations with respect to the 
fountain), and is seen near the centerline in the 
upper half of the collision zone. This is the 

point where the wall jets collide after having 
been forced up and over the separation zone at the 
base of the fountain. It is also the region of 
maximum stre amline curvature. The regions of 

maximum v'v' in the fountain are on both sides of 
the upwash and are in much the same regions as we 
expe c^ f or the shea r stress. The regions of high- 
est u u and v'v' do not overlap which suggests 
that an Important mec hanis m ma y be the inter- 
component t ransf er from u'u' to v'v". The regions 
of highest w'w' are very near the wall at the base 
of the fountain and i n the upper part of the foun- 
tain. However, w'w' is relatively large in all 
regions of turbulent flow and does not show the 
str ong spatial variations displayed by u'u' 
and v'v'. 

The turbulent shear stress normalized on the 
local centerline velocity, including the SGS con- 
tribution, Is given in Figure 13 at a few stations 
in the fountain. The SGS contribution to the 
shear stress is large in the wall jet (about 502), 
but small (less than 10%) in the fountain. The 
maximum stress is predicted to be about 0.075, 
which Is in reasonable agreement with the data of 
Saripalli, but is considerably higher than the 
value of 0.024 reported by Gilbert. The stress is 
roughly three times greater than the maximum shear 
stress of 0.022 in a "normal" planar jet [Ref. 17] 


and, therefore, consistent with the spreading rate 
which is three times greater than in the "normal" 
jet. Figure 14 gives the shear stress as a con- 
tour plot which shows regions of high stress in 
the fountain and along the edge of the separation 
bubble at the base of the turning region. 

3.3 Instantaneous Flow Field 

Flow visualization of the turbulence can 
provide insights that can never be gleaned from 
statistical data. The computer code is not cur- 
rently capable of saving all of the time dependent 
information we desire; however, we can examine an 
instantaneous three-dimensional flow field and 
learn some important points about the fountain. 

Figure 15 gives particle tracers started in 
the opposing wall jets of the instantaneous flow 
field. The tracers show a small amount of mixing, 
indicated by crossing paths, in the wall jets. In 
the fountain considerably more mixing occurs, 
especially for y > 6 which is the region of maxi- 
mum spreading rate. A mechanism that appears to 
be important is this: Blobs of fluid with higher 

u' can penetrate through the mean centerline of 
the fountain. Once through the centerline the 
blob meets less resistance to its motion since it 
is moving into a weaker flow and no adverse pres- 
sure gradient. This blob now travels on a path 
that is very different from the mean flow, which 
gives high mixing and shear stress. Another fea- 
ture of this figure is that the fountain is 
slightly tilted. The tilt might indicate a "tur- 
bulence" mechanism consisting of the fountain 
flapping back and forth. The tilt does appear in 
the mean particle paths (Fig. 8-c) which may indi- 
cate this is a problem of asymmetry induced by the 
numerics. It is also possible that a low fre- 
quency flapping motion of the fountain exists, 
which has not been removed by the time averaging; 
this would be difficult to study because of the 
long Integration times required and the dependence 
on low frequency fluctuations of the inflow condi- 
tions. 

Figure 16-a gives velocity vectors projected 
in the x-z plane in the fountain at y-12. This 
shows a region of fluid, denoted A, with a large 
negative u' component penetrating into relatively 
quiescent fluid. Mach number contours in the same 
plane, Figure 16-b, show that this is a region of 
high speed flow, indicating the vertical velocity 
is large as well (positive v'). Thus, this blob 
will contribute to negative u'v'. There is also a 
region (B) where the fountain fluid appears to be 
ingulfing a large area of ambient fluid. This may 
be the remnant of an earlier high speed blob which 
is now "mushrooming out" and forming vortices 
similar to the starting vortices from a free jet. 

The fluid interactions we have examined in 
this single "snapshot in time" Indicate that the 
fountain turbulence has a large scale structure 
caused by the penetration of blobs of high speed 
fluid through the mean fountain centerline. Data 
at other times also show these phenomena, The 
evidence for a flappii^ moti is relatively 
weak. The next step is to examine the results in 
terms that are amenable to Reynolds— averaged tur- 
bulence modeling. 
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3.4 Turbulence Mechanisms 


Evaluation of the mechanisms that contribute 
to the high turbulent stresses will be done within 
the framework of the Reynolds stress transport 
equations, which are given here in cartesian ten- 
sor notation. 
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The left hand side of Equation 4 represents the 
rate of change of the components of Reynolds 
stress along streamlines. In order, the right 
hand side terms represent: production, viscous 

dissipation, pressure strain, turbulent diffusion, 
viscous diffusion and pressure diffusion. The 
terms which are thought to be major ones will be 
given as contour plots. Identical contour levels 
and normalization are used for all terms so that 
valid comparisons of the terms can be made from 
these figures. Only resolvable turbulence contri- 
butes to these results; the SGS shear stress is 
not included. 


model the velocity triple correlations, [Ref. 18], 
which for this flow means that regions with the 
largest stress would experience the greatest loss 
in much the way heat diffuses. Comparison of 
Figures 17-c and 14 suggests that the gradient 
diffusion concept is correct near the outer edges 
of the fountain, but not in the central region of 
the fountain, between the points of maximum shear 
stress magnitude. Term IV is only slightly less 
in magnitude than the production or pressure- 
strain, and therefore relatively important. 

For the normal stresses a clear picture 
emerges regarding the roles of the different 

ter ms . Figures 18-a and 19-a show the prod uction 

of u'u' and v'v'. Production of u'u" occurs 
primarily at the collision point; this is where 
the mean flow is redirected upward, but the higher 
speed blobs penetrate through the mean center- 
l lne . The primary contributor to this terra is 
u'u' dU/dx, which is large because dU/dx is so 
large. Produc tion of v'v' occurs mainly in the 
fountain where u"v' dV/dx Is l arge. There is a 
region of negative production of v'v' at the base 
of the fountain; negative production is impossible 
to obtain with a positive definite eddy visc osity 
mode l. The pressure strain terms for u'u' 
and v'v', given in Figures 18-b and 19-b, are of 
comparable magnitude, but opposite sign, in the 
collision zone; this indicates a tr ansfe r of en- 
ergy from u'u' in the wall Jet to v'v' in the 
fountain (note that v' is streamwise in the foun- 
tain) • 

4 . SUMMARY 


The production term is important because it 
extracts energy from the mean flow and converts it 
to turbulence. The pressure-strain term cannot 
produce turbulence energy; it merely transfers 
energy among the different components. Pressure- 
strain interactions can affect the shear stress by 
changing the correlation between u' and v' . The 
objective in this examination is to determine 
which are the principal terms In the generation of 
the high shear stresses. 


The production terra for u'v' is given In 
Figure 17-a. Regions of high production are seen 
at the separation bubble at the base of the foun- 
tain and in the fountain and are nearly coincident 
with the regions of high shear stress given in 

Figure 14_. Figure 1 7— b gives the pressure-strain 

term for u'v^ and shows levels comparable to the 
production terra, but in different locations. The 
highest levels are where the wall jets first sepa- 
rate and there are moderate levels near to, but on 
either side of, the fountain centerline. The 
production and pressure-strain terms have, lo- 
cally, the same sign in most of the flow and com- 
bine to increase the level of shear stress. It 
appears that the pressure-strain term is more 
important at the base of the fountain but the 
production term is more important in the fountain. 


The gain or loss of u'v' due to turbulent 
diffusion, term IV, is given in Figure 17-c. 
Gradient diffusion models are typically used to 


Numerical simulations with the three- 
dimensional Navier-Stokes equations were used to 
study turbulence mechanisms in a VTOL upwash foun- 
tain, The primary characteristic of this flow, 
the abnormally high spreading rate of the foun- 
tain, was predicted. Large values of the Reynolds 
stresses were also predicted and these are in 
moderately good agreement with the (widely scat- 
tered) experimental data. The pressure-strain and 
production terras in the Reynolds shear stress 
transport equation have been compared and are 
shown to be of comparable magnitude and of the 
same sign. The principal mechanism for generating 
the high shear stress is the penetration of blobs 
of high speed fluid through the mean fountain 
centerline. This is reflected in the large magni- 
tudes of the prod uction terms for the late ral 
fluctuations (u'u') and the shear stress (u'v'). 

This is a flow for which advanced turbulence 
models have failed to give good results. The 
success of very large eddy simulation for what 
could be described as an engineering application 
may point the way for the prediction of other 
difficult turbulent flows. 
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Figure 1. Sketch of impinging jets with upwash 
fountain, showing computational domain. 



Figure 2. Time averaged velocity vectors in wall 
jet. 
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Figure 6* Triple velocity correlations in the 
wall jet: a) gives u'u'u'; b) gives u'u'v", c) 

gives v'v'v ' . 
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Figure 8. Mean results from unsteady fountain 
simulation: a) give velocity vectors; b) gives 

Mach number contours; c) gives streamline which 
show separation at base of fountain. 
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Figure 7, Mach number contours from a steady 
state k-e calculation, which underpredicts the 
fountain spreading rate, from Ref. 5. 
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Figure 9, Half velocity width spreading rate in 
the fountain and comparison with data from 
Gilbert, Ref, 15, 



Figure 10. Mean velocity profiles in the fountain 
and comparison with "normal" plane jet, Ref. 17. 
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Figure 12. Contour plots of normal stresses: a) 

gives u'u'; b) gives v'v'; c) gives w'w'. 


Figure 11. Streamwise normal stress, v'v', , in 
the fountain, comparison with Refs. 15 and 16; see 
Figure 10 for symbols. 
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Figure 13* Shear stress in the fountain. 



Figure 16. Instantaneous flow field in the foun- 
tain on a plane normal to the strearawise direc- 
tion: a) gives velocity vectors; b) gives Mach 

number contours. 
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Figure 15. Tracer particles in an instantaneous 
flow field (see text). 
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